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O ■ Abstract 

a • 

We present exact calculations of the zero-temperature partition function of the g-state 

^ ■ Potts antiferromagnet (equivalently the chromatic polynomial) for Mobius strips, with width 

Lj, = 2 or 3, of regular lattices and homeomorphic expansions thereof. These are compared 

^ ' with the corresponding partition functions for strip graphs with (untwisted) periodic longi- 

Y ' tudinal boundary conditions. 
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The chromatic polynomial P{G, q) counts the number of ways that one can color a graph 
G with q colors such that no two adjacent vertices have the same color (for reviews, see 
0). The minimum number of colors needed for this coloring of G is the chromatic number, 
x{G). Besides its intrinsic mathematical interest, the chromatic polynomial has an important 
connection with statistical mechanics since it is the zero-temperature partition function of 
the g-state Potts antiferromagnet (AF) |0, ^ on G: P{G, q) = Z{G, q,T = 0)paf- The Potts 
AF exhibits nonzero ground-state entropy 5*0 7^ (without frustration) for sufficiently large 
g on a given lattice graph and is thus an exception to the third law of thermodynamics (e.g. 
P). This is equivalent to a ground state degeneracy per site W > 1, since Sq = ksliaW. 
Denoting the number of vertices of G as n = v{G) and {G} = lim^^oo G, we haveQ 



W{{G},q)= hm P{G,q) 



l/n 



Let us consider strips of regular lattices with variable length and fixed width (with the lon- 
gitudinal and transverse directions taken to be the x and y directions, respectively.) An 
obvious and important question concerns the effect of boundary conditions (BC's) on the 
T = partition function. We study this effect here via new calculations of exact T = 
partition functions for arbitrarily long Mobius strips of regular lattice graphs and certain 
homeomorphic expansions thereof, together with a comparison with the corresponding par- 
tition functions for cyclic strips, i.e., strips with (untwisted) periodic longitudinal boundary 
conditions. These strips all have free transverse boundary conditions (for a comparison of 
the effects of free versus periodic transverse boundary conditions, see |jl9[)- The different 
longitudinal boundary conditions also entail different topology: the cyclic strips (topologi- 
cally ~ cylinder) are oriented, while the Mobius strips are unoriented. We use the symbols 
FBC^:, PBCx = eye, and TPBC^; = Mb. for the free, periodic, cyclic, and Mobius longi- 
tudinal boundary conditions, respectively. Some related works on calculations of chromatic 



polynomials for families of graphs are Refs. 0-||28|; see also the list in ||29|| . 

Since P{G, q) is a polynomial, one can generalize q from to C. The zeros of P{G, q) in 
the complex q plane are called chromatic zeros. Their accumulation set in the limit n — > 00, 
denoted is the continuous locus of points where W{{G},q) is nonanalytic |13|, [1^, ||, 
1^-|25]. The maximal region in the complex q plane to which one can analytically continue 



the function lydG}, q) from physical values where there is nonzero ground state entropy is 

denoted -Ri. The maximal value of q where B crosses or intersects with the (positive) real 

axis is denoted qc{{G}). 

^ At certain special points qs (typically qs = 0, 1, .., x(G)), one has the noncommutativity of limits 
\\m.q^q^ lim„^oo P(G, g)^/" 7^ Hm„^oo linV-^g^ P{G,q)^/^'', and hence it is necessary to specify the order of 
the limits in the definition of W^({G}, qs) jo). As in Ref. 1^ and our other works, we shall use the first order 
of limits here; this has the advantage of removing certain isolated discontinuities in W . 
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We label a particular type of strip graph as Gs and the specific graph of length = m 
repeated subgraph units, e.g. columns of squares in the case of the square strip and hexagons 
in the case of the kagome strip, as {Gs)m- If one thinks of the graph as embedded on 
a rectangular strip of paper, with its longitudinal ends glued with the same or opposite 
orientation, is the length of this strip in subgraph units. (Note that, e.g., for the Mobius 
square strip, if one starts at a vertex on a horizontal boundary edge and traverses a path in 
the longitudinal direction, the pathlength required to return to the original vertex is 2Lx.) 
For m greater than a minimal valueQ, the cyclic and Mobius strips of a given {Gs)m have 
the same number of vertices n, edges (bonds) e, girth g (length of minimum closed circuit 
on Gg) and number kg of circuits of length g. One has n = tsfn where tg depends on Gg- 
Writing 



n— 1 



P{{Gs),q) = Y^{-irK.,q--= (2) 
and using the results that 0, |10| h^-j = for < j < (? — 1 (whence hn = 1 and /i„_i = e) 



and hn-{g-i) = (gli) ~ ^gi it follows that for m greater than the above-mentioned minimal 
value, these hj''s are the same for the cyclic and Mobius strips of each type Gg. For a given 
Gs, as m increases, the /i„_j's for the cyclic and Mobius strips become equal for larger j (see 
below). 

A generic form for chromatic polynomials for recursively defined families of graphs, of 
which strip graphs Gg are special cases, is 

Na 

P{Gg,q) = Y,cM){aM)r (3) 

where Cj{q) and aj{q) depend on the type of strip graph Gg but are independent of m. 
Note that some Cj and aj may be nonpolynomial algebraic functions of g, in which case 
certain cancellations occur to yield the polynomial P [^U], The simplest cyclic strip is 
the circuit graph with m vertices, C^, for which P{Cm, q) = {q — 1)™" + — 1)(— 1)"*. The 
chromatic polynomials for the cyclic and twisted cyclic ladder graphs were calculated in |Q 
and have A^^^ = 4. W and B were given for these and other families in [^] (see also ||16[|). 
Various iterative methods have been used in the past to calculate chromatic polynomials for 
recursive families of graphs [0, [|T^]-[|TB|, 0, [|T^-[^. f\ A method that we have found useful 



to compute P{{Gg)m, q) for a variety of strips with FBC^ pUj , |2T| and (T)PBC2. 

^ For example, for m < 4 (to < 3) the Ly ~ 2 cyclic (Mobius) square strips exhibit special behavior 
regarding g and kg] for m larger than these respective values, they both have .g = 4 and kg = to. 

■^The method used in ||l2|, can also be used to obtain rigorous bounds on W for 2D lattices (l2[|^, 
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makes use of a generating function 

T{Gs,q,x) 



J2 P{{Gs)m,q)x 



m—mo 



(4) 



where r(Gs, g, x) = Af{Gs, q, x)/V{Gs, q, x) with 

Af{Gs,q,x) = 



x-^ 



(5) 



j=o 



and 



(6) 



deg^V degxV 

ViGs,q,x) = l+Yl ^^■(^)^' = n (1 - ^^■(^)^) 
i=i j=i 

where mo, Aj, and Xj depend on the type of strip, and Aj and bj are polynomials in q 
(mo = 2 for the square and homeomorphically expanded square strips and mo = 1 for 
the kagome strips considered below). F is calculated by means of standard theorems for 
chromatic polynomials such as the deletion-contraction theorem. The aj and Cj in (^) can 
be determined from F using eqs. (2.17)-(2.21) of (taking into account value of mo for 
a given strip). For a specified Gg, the A's in T> that contribute to P are identical to theQ 
a/s. For a given type of strip, we shall study the effect of cyclic versus Mobius longitudinal 
boundary conditions on (i) the set of aj, (ii) the W function; (iii) the locus B, (iv) the set 
of Cj, and (v) x- Quantities (i)-(iv) depend only on Gg while (v) also depends on m. f\ 

We consider first strips of the square (sq) lattice of length = m > 2 and width Ly 
vertices. The Ly = 1,2 cases have been mentioned above. For the cyclic Ly = 3 strip, P, 



W, and B were calculated in [24|. Our new result for the Ly = 3 Mobius strip is 



P{sq{Ly = 3, Mb.U q) = {q' - 3q + 1)(-1)" - (g - l)'" + (g - 2)™ - (g - 4)™ 
-(g - l)[-(g - 2)Y + k^fir + (As,j)H + (g - 1) f(As,,8)™ + (As,,9)" + (A.,,io) 



(7) 



where 

Asg,(6,7) 



(g - 2)(g2 - 3g + 5) ± |(g2 - 5g + 7)(g' - 5g^ + llg^ - 12g + 8)} 



1/2' 



For our square lattice strips, all Aj 's contribute to P. For the kagome strip and homeomorphically 
expanded square strips considered here, there are, respectvely, three A's and one A that do not contribute to 

^It is also of interest to compare chromatic polynomials for lattice strips with cyclic versus free longitudinal 
boundary conditions; see 12^, Ip, ||3|-[|5|. 
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and Xsq,j, j = 8, 9, 10, are the roots of the cubic equation 



with 



hsq,2i = - 9g + 12 
bsq,22 = q^- lOg^ + 36g2 _ 56g + 31 



(9) 

(10) 
(11) 



Jsq,23 



-{q - l)(g^ - + 29g^ - 40g + 22) . 



(12) 



For comparison, the cychc Ly = 3 resuh is p4| 

P{sq{Ly = 3, cycU, q) = {q^ - 5q' + Qq - 1)(-1) 



+ (g2-3g + l) 



(g - I)'" + (g - 2r + (g - +iq- l)[-(9 - 2) 



\2lm 



(A.,,6)™+(A.,,7)'" +(g-l) (A.,,8)"^ + (A.,,9)'" + (A.,,io) 



(13) 



We discuss these further below. 

We next consider strips of the kagome (kg) lattice comprised of m hexagons with each 
adjacent pair sharing two triangles as adjacent polygons. For the Mobius kagome strip with 
L„ = 2 we calculate 



P{kg{Ly = 2, Mb.)m, q) = -{q - 4)"^ - (g _ 1) (g - l)(g - 2 



+ 



iXk,,ir + iw^r + (9 - 1) {^k,,,r + (a..,4) 



(14) 



where 



Afc9,{l,2) 



- 6g^ + 14g2 - 16g + 10 



±<^ g» - 12g' + 64g° - 200g^ + 404g^ - 548g^ + 500g^ - 292g + 92 



1/2- 



A 



fcg,(3,4) 



q^ - + 19g - 20 ± <{ - 14g^ + 83g^ - 278g^ + 569g2 - 680g + 368 



l/2n 



(15) 



. (16) 



In contrast, the chromatic polynomial of the regular cyclic Ly = 2 kagome strip is 



P{kg{Ly = 2, cyc.U q) = (g' - 3g + l)(g - 4)™ + (g - l) (g - l)(g - 2) 
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+ 



{\kg,ir + {W2r {^^r + {^a) 



(17) 



It is also of interest to study homeomorphic expansions of the cyclic square strip.|^ We 
start with the cyclic or Mobius Ly = 2 square strip and add k — 2 vertices to each upper 
bond (edge). Define 

P{C„q) ^y ^_^y(k-V 



q{q 



s=0 



k-2-s 



q 



where P{Ck, q) = {q — l)^ + {q — 1)(^1)^ is the chromatic polynomial for the circuit graph 
with k vertices. For the twisted cyclic strip we calculate 



-P((C/l)(fc,2),Affe, 

where 

-^^,(4,5) = ; 



1)'=™ + (Dfc+2) 



Dk+1 + 



(Am)™ - {\^r (19) 



Z^fc+i - (-l)^g - 2) ± - (-l)^(g - 2)]^ + AD,i-l)'] 



1/2 



(20) 



h means homeomorphic, and {ki,k2) refers to the respective number of vertices (counting 
the end ones) on the upper and lower horizontal bounds between each vertical bond. 



For the cyclic case p4| , 
PiiCh)^k,2),cyc.,m, q) = (?' " 3g + 1)(-1)'=™ + {Dk^^r + (9 - 1) {KaT + {K^r 



(21) 



For the symmetric homeomorphic expansion, P{{Ch)(^f:,k),cyc.,m,q) and P{{Ch)(^k,k),Mb.,m,q) 
and the resultant W and B were given in p3|. Results for P, W, and B for the cyclic and 
Mobius Ly = 2 triangular strips were given in |2^ and P for the cyclic Ly = 2 triangular 
strips also in [27|. Note that for the Ly = 2 Mobius triangular strip and homeomorphically 



expanded square strip in eq. (|T^ two of the c^'s are nonpolynomial algebraic functions of q. 

We comment on how the nonpolynomial algebraic roots in the various P's yield polyno- 
mials in q. As is evident from (|^), for a given strip, the A/s arise as roots of the equation 
T> = 0. In general, T> contains some number of factors of linear, quadratic, cubic, etc. order 
in X. Consider a generic factor in V of r'th degree in x: {1 + fix + f2x'^ + . . . + frX^) , where the 
/j's are polynomials in x. This yields r A^'s as roots of the equation C,^ + + ... + = 0. 

The expressions in P involving these roots are symmetric polynomial functions of the roots, 
namely terms of the form 

r 

Sm = 5^(A^)"^ (22) 



^ A homeomorphic expansion of a graph consists in the addition of degree-2 vertices to bonds of that 
graph. 
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and, for certain Mobius strips (as well as strips with FBCx), terms of the form 



Aa — Afc 



m—l 



J^i^ar^'^'i^b)' . (23) 
fc=0 

We can then apply a standard theorem (e.g. [^) which states that a symmetric polynomial 
function of the roots of an algebraic equation is a polynomial in the coefficients, here fe, 
i = l,--,r; hence this function is a polynomial in q. For example, for Sr, one has the 
well-known formulas (due to Newton) si = — /i, S2 = /i — 2/2, S3 = —ff + 8/1/2 — 8/3, etc. 

Concerning question (i), in general we find that for a given type of lattice strip (or 
homeomorphic expansion thereof) , the terms aj in (^ are the same for the cyclic and Mobius 
strips. This can understood as follows. Recall that the linear equation for the elongation 
of the strip, together with corresponding equations for related strips produced via deletion- 
contraction operations, form the matrix equation Tv^ = Vq (eq. (2.29) of Ref. |^), where 
is the vector comprised by the generating function F for the strip of interest, together with 
generating functions for the above-mentioned related strips, and vo is a vector of chromatic 
polynomials for beginning forms of these strips. The denominator V of the generating 
function F is then given as P = det{T). The set of generating functions in the vector does 
depend on the longitudinal boundary conditions, as is clear from the fact that it has larger 
dimension for periodic than free BC^'s. However, it is the same for cyclic and Mobius BC^j's 
(the elongation is a local operation); the difference here is in the chromatic polynomials in 
Vo- For example, in the simplest case of the Ly = 2 square strip, the minimal subgraphs, viz., 
for m = 2, are a single square, or equivalently, circuit graph C4 for PBC^; but the complete 
graph K4 for TPBC^ (Mobius). [] Hence, V is the same, while A/" is different, for a given 
strip with cyclic and Mobius longitudinal BC's. Since the A/s in T> that contribute to P 
are identical to the a^-'s in eq. (^ (cf. eq. (2.20) in [^), this enables one to understand 
that these Qjj S cirG the same for the cyclic and Mobius strips of a particular type, given the 
observed property that in the cases where certain A/s do not contribute to P, these are the 
same for the cyclic and Mobius cases. 

The property that the a/s are the same for these two cases also answers questions (ii) 
and (iii). The term aj that has maximal magnitude in a region of the q plane determines 
W in that region, and hence the locus B is the solution of the equation \aj\ = \aj>\ of terms 
of maximal magnitude. It follows that W and B are the same for the cyclic and Mobius 
versions of a particular kind of strip. The definition ([^) and the equality of the W functions 
for the cyclic and Mobius strips of a given type Gg imply that, as noted above, as n increases, 
the hn-/s for these two strips become equal for progressively larger values of j. 
^The complete graph Kp is the graph consisting of p vertices all adjacent to each other. 
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On question (iv), since the Cj are functions of both the Aj's in T> and the coefficient 
functions Aj in A/", and since the latter differ for cychc and Mobius BC^-'s, the c^'s will also, 
in general, differ. However, a subset may be the same. In particular, the coefficient Cj of the 
term which is dominant in region Ri is always unity, as a consequence of the property that 
/iri = 1 in eq. @). Hence, a fortiori, this coefficient is the same for cyclic and Mobius BC^^'s. 
In cases where the leading a.j is a nonpolynomial algebraic function of g, the requirement 
that P be a polynomial implies that certain other coefficients must also be unity. Thus, for 
example, Cq and Cy are both equal to 1 for the Ly = 3 square strip, and the {\kg,i)"^ and 
{^kg,2)"^ terms both have coefficient 1 for the kagome strip. 

The situation with the chromatic number, question (v), depends on the particular strip. 
In the following, we take = m greater than the above-mentioned minimal values where 
the strips degenerate. The cyclic square strips have x = 2 for even and x = 3 for L^. 
odd, independent of Ly. For Mobius square strips, x = 2 for (L^jLy) = (e, o) or (o, e), and 
X = 3 for {Lx, Ly) = (e, e) or (o, o), where e and o denote even and odd. Both the cyclic and 
Mobius kagome strips above have x = 3, which shows that the twisting does not necessarily 
change x- The cyclic and Mobius homeomorphically expanded square strips [eqs. ( plj ) and 
(|T9|)] both have x = 3 for odd k, independent of m, while for even k, the cyclic strip has 
X = 2, 3 for m even, odd, and the Mobius strip has x = 2, 3 for m odd, even, respectively. 

For a given type of strip Gg, define C = J2f=i^j- ^'^^ ^^e Ly = 2 square strips, C = 
g(g— 1) and for the cyclic and Mobius cases H, |2^. This is also true for the Ly = 2 cyclic and 
Mobius triangular strips, respectively, as also for the cyclic and twisted homeomorphically 
expanded Ly = 2 square strips. For the Ly = 3 cyclic square strip, eq. (|l^) yields C = 
q{q — 1)^, in accord with the generalization C = P{Tr=Ly,q) = q{q — for the cyclic 

square strip of width Ly. For the Ly = 3 square Mobius strip, eq. (|^) gives C = q{q — 1). 
The Ly = 2 cyclic and Mobius kagome strips have C = q^ and C = q, respectively. 

Our comparative study thus elucidates the effect of periodic versus twisted periodic lon- 
gitudinal boundary conditions on the T = Potts antiferromagnet T = partition function 
(chromatic polynomial) for families of lattice strip graphs and homeomorphic expansions 
thereof. In the future it would also be of interest to extend this study to consider other 
combinations of longitudinal and transverse boundary conditions. 

This research was supported in part by the NSF grant PHY-97-22101. 
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